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METHODS AND APPARATUS FOR DETERMINING CARDIAC OUTPUT 

Cross-Reference to Related Application 
[0001] This application claims priority to U.S. Provisional Patent Application 
5 60/446,385, filed February 10, 2003, which is hereby incorporated by reference. 

Government Support 

[0002] This invention was made with Government Support under Contract/Grant 
Numbers 6579200 and 6890082 awarded by the National Aeronautics and Space 
10 Administration. The Government has certain rights in the invention. 

Background of the Invention 

[0003] Cardiac output (CO) is defined to be the volume of blood ejected by the heart per 
a unit time. Since CO also represents the total flow of blood supplying all the tissue beds of 

1 5 the body, it is perhaps the most indicative quantity of the state of the heart and circulation. 
CO is routinely measured in intensive care units and surgical suites in order to monitor and 
guide therapy for critically ill patients. These patients include, for example, those in shock 
(e.g., cardiogenic, hemorrhagic, or septic) or heart failure and those during and after surgery 
(e.g., coronary artery bypass grafting or heart valve replacement.) 

20 [0004] An ideal CO measurement technique would be simple to perform, inexpensive, 
noninvasive or minimally invasive and safe, and very accurate. However, none of the 
conventional measurement techniques known in the art possess all of these characteristics 
[10]. For example, the thermodilution technique, which is currently employed in most 
intensive care units and surgical suites, involves injecting cold saline into the right atrium 

25 and measuring the temperature downstream in the pulmonary artery. The average CO over 
the measurement period may than be computed from conservation of mass laws. Although 
the technique is relatively simple and inexpensive, it requires an invasive right heart 
catheterization whose safety is questionable [8,38] and is not very accurate due to the many 
assumptions upon which it is based (e.g., no saline recirculation and thorough blood mixing) 

30 [10,27]. The most accurate, conventional technique for measuring CO involves surgically 
implanting a flow probe, either electromagnetic or ultrasonic, directly on the aorta. 
Although this technique also provides a continuous measurement of CO, it requires a high 
risk thoracotomy which is rarely performed in practice. Moreover, the accuracy of the aortic 
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flow probe is highly dependent on vessel preparation and may only be accurate to within 
about 15-20 percent [10,27]. 

[0005] Although the development of an ideal CO measurement technique has proven to 
be difficult, several ideal, or near ideal, techniques are currently available for the continuous 
5 measurement of peripheral arterial blood pressure such as Finapres technology [23] and 
arterial tonometry [25]. Previous investigators have therefore sought techniques to monitor 
CO from peripheral arterial blood pressure signals. The most popular techniques in the art 
are the so-called pulse contour methods that assume the arterial tree to be well represented 
by a parallel combination of a capacitor and resistor thereby accounting for the compliance 

10 of the large arteries (AC) and the total peripheral resistance (TPR) of the small arteries. If 
the instantaneous CO supplied by the heart is represented as a current source, then the 
simple model of the heart and arterial tree in Fig. 1 A results. Most of these types of pulse 
contour methods are specifically based on mathematical formulas which are derived by 
making simplifying assumptions and approximations to this model (see, for example, [7, 19- 

15 21, 26, 40, 43, 44]). These methods have generally failed to yield good correlation between 
CO determined from analysis of ABP signals and directly measured CO over a wide range 
of physiologic conditions [39, 42]. 

[0006] Bourgeois et al [3,4] did successfully demonstrate that their pulse contour 
method when applied to an ABP signal measured centrally in the aorta could yield a quantity 

20 which varied linearly with directly measured CO (electromagnetic aortic flow probe) over a 
wide range of physiologic conditions. Their method, which makes no simplifying 
assumptions or approximations to the model of Fig. 1 A, may be explained as follows. 
During the diastolic period of each cardiac cycle, the heart is filling and not supplying blood 
to the arterial tree (see Fig. IB). Thus, according to the model, ABP decays with a time 

25 constant t d equal to the product of TPR and AC during each diastolic period (see Fig. IB). 
Since AC is essentially constant over a wide pressure range and on the time scale of days [4, 
18,37], CO could then be monitored to within a constant scale factor (equal to the reciprocal 
of AC) by dividing the ABP by x D . 

[0007] Bourgeois et al specifically demonstrated that beat-to-beat CO may be 
30 monitored from x D and the governing differential equation of the model of Fig. 1 A. Thus, a 
key step of the pulse contour method of Bourgeois et al is to fit an exponential to the 
diastolic decay portion of an ABP wavelet in order to measure x D . Osborn et al [34] 
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introduced essentially the same method prior to Bourgeois et al but their experimental 
validation was not as complete or compelling. 

[0008] Bourgeois et al were able to validate their pulse contour method with respect to 
a canine ABP signal measured centrally in the aorta, because the diastolic portion of such a 
5 signal usually resembles an exponential decay (see Fig. 2A). These investigators 

specifically identified the position in the aorta at the level of the dorsal insertion of the 
diaphragm as the optimal site for observing an exponential diastolic decay. However, 
Bourgeois et al acknowledged that central ABP is rarely obtained clinically because of the 
difficulty in inserting a catheter retrogradely via a peripheral arterial blood vessel and the 

10 risk of blood clot formation and embolization. Moreover, they recognized that, in peripheral 
ABP signals which are routinely made available in intensive care units and surgical suites 
usually via a more simple and safe radial artery catheterization, an exponential diastolic 
decay is usually not apparent (see Fig. 2B). The method of Bourgeois et al therefore cannot 
generally be applied to readily available peripheral ABP signals. In fact, its application to 

15 central ABP signals may be somewhat limited, as Cundick et al [9] reported that they could 
not identify an optimal location in the human aorta in which the diastolic portion of the ABP 
signal appeared as a pure exponential decay. 

[0009] Other pulse contour methods that are based on more complex representations of 
the arterial tree have also been developed (see, for example, [5, 1 1, 12, 16, 24, 31-33, 50]). 
20 However, these techniques required the analysis of one, or even two, central ABP signals. 
Thus, their clinical utility is also severely limited. 

[0010] Several techniques have more recently been introduced in an attempt to monitor 
CO from ABP signals measured peripherally. Techniques based on an adaptive aorta model 
which require ABP signals measured at two peripheral sites - the carotid artery and the 

25 femoral artery - have been developed [36, 46]. However, catheters are usually not placed 
for prolonged periods of time at either of these sites in intensive care units or surgical suites 
due to issues of safety. Another previous technique is based on an empirically-derived 
formula which involves the calculation of the derivative of the ABP signal [14]. However, 
in order to mitigate the corruptive effects of wave reflections on the derivative calculation, 

30 this technique also requires two peripheral ABP measurements, one of which is obtained 
from the femoral artery. Other techniques based on a learning approach have been 
previously proposed [6, 15, 30]. However, these techniques require extremely large sets of 
training data consisting of simultaneous measurements of CO and ABP signals obtained 
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over the entire range of physiologic conditions. Moreover, the success of these techniques 
was only demonstrated with central ABP signals or only over a narrow physiologic range. 
Finally, Wesseling et al [1, 48, 49] and Linton and Linton [28] have recently proposed 
model-based techniques which require only the analysis of a single radial artery pressure 
5 signal. However, Linton and Linton showed that their technique was reasonably accurate 
only over a narrow range of physiologic conditions, and several previous studies have 
demonstrated the inadequacy of the method of Wesseling et al (see, for example, [13, 22]). 
[0011] It is evident that there remains a need in the art for methods and apparatus for 
determining CO reliably and accurately using information obtained from the arterial blood 
10 pressure signal. In particular, there remains a need in the art for methods and apparatus for 
determining CO reliably and accurately using information obtained from the peripheral 
arterial blood pressure signal. 

Summary of the Invention 

1 5 [0012] The present invention addresses this need, among others. In one aspect, the 
invention provides for the measuring of a physiologic signal indicative of cardiovascular 
system activity, e.g., an arterial blood pressure (ABP) signal over a plurality of cardiac 
cycles. From analysis of the arterial blood pressure signal or other physiological signal the 
times of the cardiac contractions are identified. Then the relationship between the times of 

20 the cardiac contractions and the physiologic signal over a plurality of cardiac cycles is 
mathematically analyzed. From this analysis a mathematical relationship between the 
occurrence of the cardiac contraction at a certain time and the subsequent time evolution of 
the arterial blood pressure over a time period greater than one cardiac cycle is obtained. 
This mathematical relationship is then used to determine a dynamical property of the 

25 system. For example, in one preferred embodiment the mathematical relationship is the 
impulse response function between the occurrence of cardiac contractions and the arterial 
blood pressure signal. In one preferred embodiment the dynamical property is the time 
constant which describes the decay of the impulse response function over a time interval, 
e.g., a time interval such as between 2 and 4 seconds following the maximum height of the 

30 exponential. In another preferred embodiment the dynamical property is the impulse 
response function itself. 

[0013] In particular, in a preferred embodiment the invention provides a method and 
apparatus for monitoring CO by analyzing the long time scale variations (greater than a 
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cardiac cycle) in a single ABP signal, which may be obtained at any site in the systemic or 
pulmonary arterial tree. The present invention determines t d through the analysis of long 
time intervals (60-90 second intervals in a preferred embodiment) of an ABP signal 
according to the following three steps (see Fig. 3). The first step involves constructing an 

5 impulse train signal, x(t), in which each impulse is located at the start of a cardiac 
contraction and has an area equal to the ensuing arterial pulse pressure (or an arbitrary 
constant value in another preferred embodiment). The constructed signal therefore 
approximately reflects cardiac contractions in terms of timing and output (or just timing, in 
those embodiments in which the area of the impulses is an arbitrary constant). 

1 0 [0014] The second step deals with determining the dynamical properties of the arterial 
tree through the characterization of the relationship between the cardiac contractions and the 
ABP signal. This is achieved by estimating an impulse response function (h(t)) which when 
convolved with x(t) best fits the ABP signal (y(t)\ according to any of a number of criteria. 
The estimated h(t) represents a normalized ABP response to a single cardiac contraction. 

15 The final step involves fitting an exponential to the tail end of the diastolic decay portion of 
the estimated h(t) in which the faster wave reflections have vanished in order to determine 
t D . Accurate determination of t d is achieved by virtue of h(t) coupling the long time scale 
variations 'mx(t) toy(t). 

[0015] The present invention, which in a preferred embodiment includes an analog-to- 
20 digital converter, a buffer system, a signal processing unit, and a display and alarm system 
(see Fig. 4), may thus be utilized to monitor CO and/or TPR, to within constant scale factors 
from any measured ABP signal despite the presence of wave reflections. Note that absolute 
CO and/or TPR of the systemic or pulmonary arterial tree (depending on the signal 
measurement site) may also be determined, if desired, by calibration with a single, absolute 
25 measure of CO such as a thermodilution measurement. 

[0016] The present invention may be employed in intensive care units and surgical 
suites in which invasive radial ABP signals are routinely available and pulmonary ABP 
signals are sometimes available. The present invention may also be applied to 
noninvasively measured peripheral ABP signals (e.g., Finapres technology, arterial 
30 tonometry) or noninvasively measured signals related to ABP (e.g. , fingertip 

photoplethysmography [45], ear densitography [17]). Thus, the present invention could 
easily be employed in primary care settings, emergency rooms, and regular hospital beds in 
order to facilitate the evaluation of the patient's heart and circulatory state. 
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[0017] This application refers to various patents and publications. The contents of all of 
these are incorporated by reference. 

Brief Description of the Drawing 

5 [0018] Figure 1 A is a simple, lumped model of the heart and arterial tree (either 
systemic or pulmonary). CO is cardiac output; ABP, arterial blood pressure; TPR, total 
peripheral resistance; and AC, arterial compliance. 

[0019] Figure IB is a plot of a realistic, continuous CO signal and a plot of the ABP 
signal response to the CO signal determined according to the simple model of Figure 1 A. 
10 The time constant x D governs the dynamical behavior of the simple model. 

[0020] Figure 2 A is an example of a central ABP signal measured in an experimental 
sheep preparation. 

[0021] Figure 2B is an example of a peripheral ABP signal measured in the brachial 
artery of an experimental dog preparation. 
1 5 [0022] Figure 3 is an illustration of how one embodiment of the present invention 

determines x D from a peripheral ABP signal. PP is arterial pulse pressure; R, the time of the 
R-wave in a surface electrocardiogram; x(t), an impulse train signal representing cardiac 
contractions; y(t), an ABP signal; and h{t\ the impulse response which when convolved with 
jc(0 best fits y(t). 

20 [0023] Figure 4 is a block diagram of one embodiment of the present invention. 

[0024] Figure 5 is the result of a pilot experimental rabbit study which demonstrates the 
correspondence between relative CO changes determined from a femoral ABP signal 
according to the present invention and an aortic flow probe. CO was altered by various 
interventions (e.g., pacing, inferior vena cava balloon occlusion, and nitroglycerine). The 

25 solid line is the line of identity, and the dashed line is the line that best fits the data points. 
[0025] Figure 6 shows example excerpts of measured (a) central, (b) "radial", and (c) 
femoral arterial blood pressure (ABP) signals from an experiment performed in swine. 
[0026] Figure 7 presents results from animal 1 in an experiment in which the method of 
the invention was tested in swine. The top panel illustrates the aortic flow probe cardiac 

30 output (CO) trend (solid), estimated and calibrated CO trends from the "radial" (dash) and 
femoral (dash-dot) arterial blood pressure (ABP) signals, and the intervention and duration 
(underline). The second panel depicts x D trends estimated from the "radial" (dash) and 
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femoral (dash-dot) ABP signals. The bottom two panels illustrate the mean ABP (MAP) 
and heart rate (HR) trends. 

[0027] Figure 8 present results from animal 5 in an experiment in which the method and 
apparatus of the invention were tested in swine. Panels correspond to those in Figure 7. 

5 

Detailed Description of Certain Preferred Embodiments of the Invention 

[0028] The present invention encompasses the recognition that there is significant 
information present in the ABP signalwhen measured over long time scales (greater than a 
single cardiac cycle), and that this information, referred to herein as "long time scale 
10 information", may be used to facilitate determination of dynamical properties of the 
systemic or pulmonary arterial tree. 

[0029] Long time scale information incorporates information reflecting the variability of 
the ABP signal between beats in addition to, or instead of, reflecting only information about 
the ABP signal within single beats. Although long time scale information is acquired by 

1 5 measuring the ABP signal over a period greater than a single cardiac cycle (typically 
seconds to minutes), not all quantities derived from measurements performed over time 
periods greater than a single cardiac cycle include long time scale information. For 
example, it is common to measure ABP over a plurality of cardiac cycles and average the 
measured amplitudes to obtain mean (average) ABP. However, mean ABP does not 

20 incorporate long time scale information, because it contains no information reflecting the 
variability between beats. Once an average value is found (e.g., by integrating the ABP 
signal over a time interval and dividing by the length of the time interval), information 
reflecting differences in length and/or amplitude of the ABP signal is lost. Thus the mean 
ABP does not incorporate long time scale information. 

25 [0030] Capturing the information present over long time scales offers a number of 
advantages for the determination of dynamical properties of the systemic or pulmonary 
arterial tree. In particular, utilizing such information provides a means of accurately 
determining such properties through measurement of the ABP signal, despite the fact that 
the ABP signal (particularly the peripheral ABP signal) is corrupted by wave reflections that 

30 occur at sites of impedance mismatch (e.g., vessel bifurcations). 

[0031] The inventors have recognized that at sufficiently long time scales, the 
wavelengths of the propagating waves may be longer than the dimension of the arterial tree. 
Thus the wave reflections corrupt peripheral ABP signals, but only on relatively short time 
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scales which, however, may exceed the duration of a single cardiac cycle, leaving the signal 
on longer time scales relatively undisturbed. This implies that the lumped model in Fig. 1 A 
is a valid representation of the long time scale dynamics of the arterial tree (either systemic 
or pulmonary) despite its limitations as an accurate representation of the arterial tree over 
5 single cardiac cycles. Moreover, the arterial tree is continuously being excited at time scales 
greater than a cardiac cycle by ongoing perturbations (eg., breathing) and the dynamic, 
compensatory response of the regulatory system (e.g., arterial baroreflex) [2]. Thus 
significant long time scale information is generally present in ABP signals, and the present 
invention utilizes this information to provide methods for determining dynamical properties 

10 of the systemic or pulmonary arterial tree. 

[0032] The invention provides a method for determining a dynamical property of the 
systemic or arterial tree comprising steps of: (a) measuring a physiologic signal over a 
plurality of cardiac cycles; (b) obtaining a relationship between the timing of a cardiac 
contraction and the evolution of the physiologic signal over a time period greater than that of 

1 5 a single cardiac contraction cycle by analyzing the physiologic signal over a plurality of 
cardiac cycles; and (c) using the relationship to determine the dynamical property. In 
general, the physiologic signal is a signal indicative of cardiovascular system activity. For 
example, in certain embodiments of the invention the physiologic signal is an arterial blood 
pressure (ABP) signal. In certain embodiments of the invention the physiologic signal is a 

20 signal related to the ABP signal. Such related signals include, but are not limited to, 

pressure signals obtained using fingertip photoplethysmography [45], ear densitography 
[17], etc. In certain embodiments of the invention the physiologic signal is the arterial - 
systemic filling pressure difference (ASFPD). For purposes of description, the invention 
will be described in terms of the ABP signal, but it is to be understood that other physiologic 

25 signals may be similarly used. 

[0033] In general, a dynamical property of a system is a characteristic of the system that 
relates to how the system responds over time to a change in one or more of the parameters of 
the system. For example, a mathematical expression that relates future values of one or 
more signals generated by a system to past and present values of signals either sensed or 

30 generated by the system would constitute a dynamical property of the system. A 

characteristic time is also a dynamical property. In general, a characteristic time indicates 
the time scale of the temporal evolution of a function. (For example, for an exponential 
function given by y(t) = the time constant t is a characteristic time of the function.) 
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[0034] In particular an impulse response function which enables one to compute the 
expected future values of a signal generated by the system from past values of that signal or 
past values of other signals either generated or sensed by the system would constitute a 
dynamical property of the system. The characteristic time of the decay of that impulse 
5 response function would also be a dynamical property of the system. Note that the 

dynamical property need not fully predict how the system responds over time to a change in 
one or more of the parameters of the system, rather it need only be descriptive of the 
response. Dynamical properties of systems can often be estimated from analysis of signals 
associated with the system. 

10 [0035] The relationship between the timing of a cardiac contraction and the evolution of 
the ABP signal over time may be obtained in a number of different ways. Determination of 
one or more dynamical properties of the systemic or pulmonary arterial tree in turn allows 
one to obtain values for a variety of important parameters that characterize the 
cardiovascular system, including, but not limited to, cardiac output, total peripheral 

15 resistance, cardiac index, stroke volume, characteristic time constant, etc. 

[0036] One application of the methods described above is in measuring cardiac output 
(CO). The invention provides a method of determining cardiac output to within a constant 
scale factor comprising steps of: (a) measuring a physiologic signal over a plurality of 
cardiac contraction cycles; (b) estimating a function that represents the response of the 

20 physiologic signal to a cardiac contraction over a time period greater than that of a single 
cardiac cycle; (c) determining a characteristic time of the function; (d) determining cardiac 
output to within a constant scale factor by dividing the magnitude of the physiologic signal 
by the characteristic time obtained in step (c). The invention further provides a method of 
determining total peripheral resistance to within a constant scale factor comprising steps of: 

25 (a) measuring a physiologic signal over a plurality of cardiac contraction cycles; (b) 
estimating a function that represents the response of the physiologic signal to a cardiac 
contraction over a time period greater than that of a single cardiac cycle; and (c) determining 
a characteristic time of the function, wherein total peripheral resistance is given to within a 
constant factor by the characteristic time. 

30 [0037] In certain embodiments of the invention rather than determining the 

characteristic time of the function estimated in part (b) of the above methods, a second 
function that represents the response of a different physiologic signal to a cardiac 
contraction over a time period greater than that of a single cardiac contraction is estimated, 
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and the characteristic time of this second function is determined and used in step (d). 
Accordingly, the invention provides a method of determining cardiac output to within a 
constant scale factor comprising steps of: (a) measuring a first physiologic. signal over a 
plurality of cardiac contraction cycles; (b) measuring a second physiologic signal over a 

5 plurality of cardiac contraction cycles; (c) estimating a function that represents the response 
of the second physiologic signal to a cardiac contraction over a time period greater than that 
of a single cardiac cycle; (d) determining a characteristic time of the function; and (e) 
determining cardiac output to within a constant scale factor by dividing the magnitude of the 
first physiologic signal by the characteristic time obtained in step (d). The methods will now 

10 be described in more detail. 

[0038] A feature common to the techniques discussed above for monitoring CO from 
continuous ABP is that the signal analysis is considered only within individual cardiac 
cycles. Because of the presence of wave reflections at these time scales, these techniques 
are limited in that they 1) can only be applied to central ABP signal in which the cumulative 

1 5 effects of the pulse reflections are largely attenuated; 2) necessitate two peripheral ABP 
signal measurements which are not usually obtained clinically; 3) require a large set of 
training data obtained over a wide range of physiologic conditions, or 4) are reasonably 
accurate only over a limited physiologic range. 

[0039] In particular, methods that attempt to determine CO by measuring the time 
20 constant (x D ) of the arterial tree, which could then be used to compute CO to within a 

constant scale factor (equal to the reciprocal of AC) have had only limited success for the 
following reasons. In peripheral ABP signals such as those typically available in intensive 
care units and surgical suites via a technique such as radial artery catheterization (and also 
peripheral ABP signals obtained via non-invasive techniques such as fingertip 
25 photoplethysmography or ear densitography), an exponential decay is not usually apparent 
(see Fig. 2B). This is because, as mentioned above, the arterial tree is not a lumped system 
but actually a distributed system with impedance mismatches throughout the system due to 
vessel tapering and bifurcations as well as changes in vessel caliber. Thus the diastolic 
portion (as well as the systolic portion) of peripheral ABP signals is corrupted by wave 
30 reflections that occur at each site of impedance mismatch. (Note that the complexity of 
these sites as well as their varying distances from the aorta result in reflected waves with 
large phasic differences which generally tend to mitigate the cumulative effects of these 
waves on the central ABP signal.) Even in the case of central APB signals it may not be 
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possible to identify an optimal location in the human aorta at which the diastolic portion of 
the ABP signal appears as a pure exponential decay [9]. Similarly, efforts to locate a 
suitable position in the pulmonary arterial tree in which the diastolic portions of the 
pulmonary ABP signal appear as an exponential decay have been unsuccessful [41]. 
5 [0040] Thus the presence of wave reflections has impeded development of accurate 
methods to determine the time constant of the arterial tree, which could then be employed to 
compute CO. As desribed below, the present invention overcomes this difficulty through 
the use of long time scale information to accurately estimate td, from which CO may then 
be determined using the following formula: 

10 

[0041] CO = (ABP]£AC) 

[0042] The scale factor, AC, for a particular individual may be determined by obtaining 
15 a single absolute measurement of CO, e.g., by thermodilution, and then solving for AC in 
the above formula using the estimated value for Id- Alternately, AC may be estimated using 
tables or nomograms, which are well known in the art and may be based on parameters such 
as age, weight, height, or particular disease status (see, e.g., U.S. Patent No. 6,485,431). 
Additional parameters such as total peripheral resistance (TPR), cardiac index, stroke 
20 volume (SV), etc., may also be determined using well known relationships. For example, 
TPR = (t d )/AC, TPR = ABP/CO, and SV = (CO)/heart rate. 

[0043] The CO signal estimated from the formula given above will reflect true CO over 
time periods greater than or equal to a single cardiac cycle, but will generately not 
accurately reflect the cardiac flow signal within a single cardiac cycle. For this reason, it 

25 may be desirable to average the ABP signal or the estimated CO signal over each cardiac 
cycle. Alternatively, in certain embodiments of the invention the ABP signal or estimated 
CO signal is filtered using a low-pass filter with a characteristic response time greater than 
the duration of a typical cardiac cycle. Another alternative is to simply average the ABP 
signal or estimated CO signal over a time period longer than the duration of a typical cardiac 

30 cycle (multibeat averaging). In general, a single-beat average approach or by use of a low- 
pass filter with a fairly short characteristic time would be expected to retain the most 
information regarding time variation of the CO signal as compared with simple averaging of 
the CO signal over long time periods. The time constant t d changes slowly in time 
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(because over a time scale of tens of seconds peripheral resistance is slowly varying and the 
arterial compliance AC may be regarded as essentially constant), however the cardiac output 
itself can vary much more rapidly - on a beat-to-beat basis. Thus it is generally preferable to 
estimate the CO from the formula given above together with use of a single beat average or 
5 low-pass filtering approaches - even if the constant td is estimated from long epochs of data 
comprising perhaps tens of seconds - as opposed to simple averaging over time scales long 
compared to the duration of a single cardiac cycle. 

[0044] Although absolute CO may be determined using the formula given above, one 
important aspect of the invention is the recognition that in many circumstances it is not 

10 necessary to obtain a value for absolute CO. For example, in the context of continuous 

monitoring in the acute setting (e.g., in intensive care units), it is changes in CO rather than 
absolute CO that is most clinically relevant. Thus determination of the proportionality 
constant is unnecessary, and this potential source of error may thus be avoided. In other 
words, the present invention may be used to monitor CO (e.g., identify and quantify changes 

1 5 in CO) instead of (or in addition to) determining absolute CO. 

[0045] According to the method of the invention an analog ABP signal is measured 
invasively or noninvasively at any site in the systemic or pulmonary arterial tree. The 
analog signal is quantized and sampled. For example, in a preferred embodiment of the 
method the signal is quantized at 12 bits and sampled at 90 Hz. It is noted that these values 

20 are exemplary only, and one of ordinary skill in the art will readily be able to select other 
appropriate values. A signal representing cardiac contractions is constructed through the 
formation of an impulse train in which each impulse is located at the start of a cardiac 
contraction and has an area equal to the ensuing arterial pulse pressure, i.e., the pulse 
pressure that results from that cardiac contraction (see Fig. 3). The start of each cardiac 

25 contraction is determined by detecting the onset of the upstroke of each ABP wavelet, which 
may be done by any of a variety of methods known in the art, e.g., those described in [51], 
[52], or [53]. The arterial pulse pressure for each wavelet is given by the difference in the 
maximum value of ABP and the value of ABP at the onset of the upstroke. 
[0046] Alternatively, the area of each impulse may be set to an arbitrary constant value. 

30 In another preferred embodiment, a surface electrocardiogram (ECG) is measured 

simultaneously with the ABP signal. The two signals are quantized at 12 bits and sampled 
at 360 Hz. It is noted that these values are exemplary only, and one of ordinary skill in the 
art will readily be able to select other appropriate values. The start of each cardiac 
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contraction may then be established by detecting each R-wave of the ECG. In certain 
embodiments of the invention the constructed impulse train and ABP signal are then 
decimated, e.g., for purposes of noise reduction. For example, in one embodiment of the 
invention the constructed impulse train and ABP signal are decimated from 360 Hz to 
5 90 Hz. 

[0047] The single contraction ABP response (normalized by approximately the average 
arterial pulse pressure when each impulse is scaled to the ensuing arterial pulse pressure), 
which quantitatively characterizes the dynamical properties of the arterial tree, is determined 
by estimating the discrete-time impulse response function (h(t)) which when convolved with 
10 the 90 Hz impulse train (x(t)) "best" fits the 90 Hz ABP signal (y(t)) in the least squares 
sense (see Fig. 3). 

[0048] The impulse response function is assumed to be well represented by an 
autoregressive moving average (ARMA) model which is given below: 

m n 

y(t) = aiy tt ~ *) + Yl ^ ~ *) + e (*)' 
15 i=l i=1 

where e(t) is the residual error term, m and n limit the number of terms in the model (model 
order), and the set of parameters {aj, bi\ completely defines h(t) [29]. Because the ARMA 
model is parametric, causality may be imposed, which is necessary for reliably estimating 
20 h(t) as x(t) and y(t) are related in a closed-loop fashion (i.e., y(t) also influences x(t) through 
the autonomically mediated heart rate baroreflex) [47]. 

[0049] For a fixed model order, the set of parameters is estimated from 60-90 second 
intervals of x(t) and y(t) through the least-squares minimization of the residual error term, 
which has an analytic solution [29]. The model order is determined by an ARMA parameter 

25 reduction algorithm that penalizes for the degree of model complexity [35] . Prior to 

estimation of h(t), x(t) and y(t) may be lowpass filtered in order to amplify the contribution 
of long time scale energy such that the least squares fit between x(t) and y(t), at these time 
scales, is prioritized. Note that any other parametric model (e.g., autoregressive moving 
average with exogenous input (ARMAX) model [29]) may be employed in various 

30 embodiments of the invention to represent the structure of h(t), and any other minimization 
criterion (e.g., absolute error) may be utilized to find the "best" h(t). In certain embodiments 
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of the invention the number of parameters in the model is selected at least in part based on 
the heart rate. 

[0050] The t d quantity is determined by finding the "best" exponential that fits h(t) over 
a selected time interval following the time of the maximum value of h(t), preferably a time 
5 interval in which the faster wave reflections have become minimal. For example, as shown 
in the plot of h(t) presented in Figure 3, the contribution of the faster wave reflections 
becomes minimal at approximately 1.5-2 seconds, as evidenced by the decrease in 
fluctuations of h(t). Following this time the impulse response may be accurately 
approximated as an exponential. Thus in certain preferred embodiments of the invention the 

10 selected time interval begins approximately 1.5 seconds following the time of maximum 
h(t), more preferably approximately 2 seconds following the time of maximum h(t). For 
example, it has been found that a time interval of 2 to 4 seconds following the time of 
maximum value of h(t) is suitable. Longer time intervals may also be used. Typically the 
appropriate time interval is predetermined, but in certain embodiments of the invention it 

1 5 may be selected as the measurements are being made. 

[0051] The following exponential equation is the basis of the resulting least squares 
problem where A and t d are parameters to be estimated through the least squares 
minimization of w(t). 

2 0 h(t) = Ae^ +w(t), 

[0052] By first log transforming h(t) over the interval of interest (which is always 
greater than zero), the optimal estimate of the parameters in the least squares sense may be 
estimated through an analytic linear least squares solution [4]. CO may then be computed to 
25 within a constant scale factor equal to 1/AC through the ratio of the ABP signal to id as 
discussed above. Note that TPR of the systemic or pulmonary arterial tree (depending on 
the signal measurement site) is trivially given, to within a constant scale factor equal to AC, 
by x D . 

[0053] In the embodiment of the invention described above proportional CO is 
30 determined by dividing ABP by t d . In another embodiment of the invention, rather than 

employing ABP, the method uses the arterial - systemic filling pressure difference (ASFPD). 
The ASFPD is determined by subtracting systemic venous pressure from arterial blood 
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pressure [54]. The systemic filling pressure can either be measured, or more commonly 
estimated. All the analyses described above can then be performed on the ASFPD rather 
than on the ABP signal, including estimation of the cardiac output. In this embodiment the 
ASFPD is divided by the time constant to obtain a signal proportional to cardiac output. 
5 [0054] This embodiment of the invention offers a number of potential advantages. In 
the absence of cardiac contractions the arterial blood pressure would decay over time to the 
systemic filling pressure, thus it is not fully accurate to describe the decay of the ABP 
impulse reponse function at long times as an exponential that decays to zero. However, 
since in the absence of cardiac contractions the ASFPD does decay to zero, a description in 

10 terms of an exponential decay at long times is more appropriate for the impulse response 
function of the ASFPD. Furthermore, since the time constant of the impulse response 
function is estimated at long times when the impulse response function has already decayed 
substantially, even if the systemic filling pressure is small, there may be a significant 
difference in the decay constant estimated from the impulse response functions of the ABP 

15 compared to the ASFPD - with the analysis of the ASFPD providing a more accurate result. 
[0055] Fig. 4 depicts a block diagram illustrating a preferred embodiment of the present 
invention. An analog ABP signal is fed into an analog-to-digital converter as it is being 
measured. The ABP signal may be acquired using standard methods, such as those 
mentioned above. In certain embodiments of the invention a surface electrocardiogram is 

20 also obtained (ECG), e.g., via standard ECG leads. The digitized ABP signal is stored in a 
buffer system. The most recent 60-90 second intervals of the digitized signal are transferred 
from the buffer system to a processing unit which analyzes the signal according to Fig. 3 in 
order to estimate x D . The buffer and processing unit may be implemented using, for 
example, any standard microcomputer running appropriate software to implement the 

25 mathematical operations described above. The software components of the invention may 
be coded in any suitable programming language and may be embodied in any of a ragne of 
computer-readable media including, but not limited to, floppy disks, hard disks, CDs, zip 
disks, DVD disks, etc. Outputs such as proportional CO, CO, and TPR may be displayed on 
a visual display such as a computer screen and/or may be printed or transmitted to a remote 

30 location. The ECG, and analysis thereof, may also be displayed. In a preferred embodiment 
of the system the process is continuously repeated thereby providing the on-line monitoring 
of CO and TPR (with a delay of 30-45 seconds). Finally, in certain embodiments of the 
invention an alarm is triggered upon excessive decreases in CO. 
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[0056] Alternatively, the AC proportionality constant may be computed with a single 
absolute measure of CO (e.g., thermodilution) through the product of t d and the measured 
CO divided by ABP, as discussed above. The proportionality constant AC may then be 
utilized to obtain absolute measures of CO and TPR. Note that physiologic changes in AC 
5 due to disease or aging may also be monitored with multiple, simultaneous measurements of 
absolute CO and an ABP signal. 

[0057] The present invention was evaluated in a pilot study in which femoral ABP and 
aortic flow probe CO were simultaneously measured in an experimental rabbit preparation 
during various interventions known to alter CO (e.g., pacing, inferior vena cava balloon 

10 occlusion, and nitroglycerine). Fig. 5 illustrates that the resulting relative changes in CO 
with respect to baseline predicted by the present invention correspond to those changes 
determined from the aortic flow probe. In Fig. 5 the x-axis represents the percentage change 
in CO as measured using the aortic flow probe, while the y-axis represents the percentage 
change in CO as determined using the present invention. The solid line (line of identity) 

1 5 represents results that would be achieved by a technique that precisely replicated the aortic 
flow probe measurements. The dashed line represents the values measured using the present 
invention and demonstrates a very close fit.) 

[0058] The invention was further evaluated in six experiments in swine, in which 
peripheral ABP signals and independent CO via an aortic flow probe were simultaneously 

20 measured over a wide physiologic range. Six Yorkshire swine (30-34 kg) were studied 

under a protocol approved by the MIT Committee on Animal Care. The animals were given 
intramuscular tiletamine-zolazepam, xylazine, and atropine prior to endotracheal intubation. 
The swine were then maintained in a deep plane of anesthesia with inhaled isoflorane 0.5%- 
4%. Positive-pressure mechanical ventilation at a rate of 10-15 breaths/min and a tidal 

25 volume of 500 ml was employed. 

[0059] Physiologic transducers were placed as follows. 7.5 French sheath introducers 
(Arrow International, Reading, PA) were placed in the bilateral femoral arteries. A 
micromanometer-tipped catheter (SPC 350, Millar Instruments, Houston, TX) was fed 
retrograde to the thoracic aorta from the femoral artery for central ABP. The catheter was 

30 specifically positioned to achieve a diastolic decay that appeared as exponential as possible. 
The second introducer was attached to stiff fluid-filled tubing (Arrow International) and an 
external pressure transducer (TSD104A, Biopac Systems, Santa Barbara, CA) for femoral 
ABP. The chest was opened with a midline sternotomy. An ultrasonic flow probe was 
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placed around the aortic root for gold standard CO (T206 with A-series probes, Transonic 
Systems, Ithaca, NY). Finally, a 23- or 25-gauge angiocatheter was placed as distal as 
possible to the brachial artery and attached to an external pressure transducer via short, rigid 
tubing for "radial" ABP. Each transducer output was interfaced to a microcomputer via an 
5 A/D conversion system (MP150WSW, Biopac Systems). The data were recorded at a 
sampling rate of 250 Hz and 16-bit resolution. 

[0060] In each animal, a subset of the following interventions was performed over the 
course of 75 to 150 minutes to vary CO and other hemodynamic parameters: infusions of 
volume, phenylephrine, dobutamine, isuprel, esmolol, nitroglycerine, and progressive 

10 hemorrhage. Several infusion rates were implemented followed by brief recovery periods. 
[0061] The technique was applied off-line to six-minute intervals (overlapping by three 
minutes) of the digitized "radial" and femoral ABP signals to estimate t and proportional 
CO trends. The corresponding gold standard CO trends were established by averaging the 
aortic flow probe measurements over the identical time intervals. Gold standard x trends 

1 5 were similarly sought by applying the technique of Bourgeois et al [4] to the central ABP 
signals. As a metric for comparison between an estimated trend ( X(i) ) and the 
corresponding gold standard trend (^ 0 (0 ), the root-mean-square-normalized-error 
(RMSNE) in percent was computed as follows: 



20 where N represents the number of analyzed six-minute intervals and the argument i denotes 
the i th analyzed six-minute interval. In order to use this metric to compare an estimated 
proportional CO trend with the corresponding absolute gold standard CO trend, the former 
trend was first scaled to have the same mean value as the latter trend for each animal. The 
correlation coefficient (p) between the estimated and gold standard trends was also 

25 calculated as a scale-invariant metric for comparison. 

[0062] Figure 6 illustrates example excerpts of the digitized central, "radial", and 
femoral ABP signals. The central ABP excerpt is from the beginning of the recording 
period in which the aortic catheter was positioned such that the diastolic decay appeared as 
exponential as possible. Unfortunately, the diastolic decay of the central ABP signal did not 

30 consistently appear as a single exponential throughout the recording period in which various 
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interventions were employed. The "radial" and femoral ABP excerpts are from the same 
time but later in the recording period. Although the diastolic decay of the femoral ABP 
excerpt appears smooth, it cannot be adequately represented by a single exponential 
function. 

5 [0063] Table 1 summarizes the results for each animal, and Figures 7 and 8 are 

examples of the corresponding trends for animal 1 (worst result) and animal 5 (best result). 
(The femoral ABP signal of animal 4 was not analyzed due to excessive damping.) The 
table and figures illustrate the wide physiologic range imposed by the interventions and 
strong agreement between the estimated and gold standard CO trends (in terms of CO 

1 0 RMSNE and visually). This strong agreement was confirmed by high overall r values 
(mean±SD) between the gold standard CO trends and the CO trends estimated from the 
"radial" (0.84±0.07) and femoral (0.86±0.05) ABP signals. Additionally, Table 2 shows 
that the CO errors (difference between calibrated, estimated CO trends and corresponding 
gold standard CO trends) were largely uncorrected with CO, mean ABP (MAP), and heart 

15 rate(HR). 
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[0064] 

Table 1. Summary of results. CO is cardiac output; MAP, mean arterial pressure; HR, heart 
rate; RMSNE, root-mean-square-normalized-error; and t d , long time constant of arterial 
tree. 

5 



ANIMAL 


CO RANGE 
[L/MIN] 


MAP RANGE 
[MMHG] 


HR RANGE 
[BPM] 


CO RMSNE [%] 
femoral "radial" 


t d RMSNE 
[%] 


1 


1.6-5.2 


29-100 


96-180 


19.9 


19.1 


8.3 


2 


2.3-4.2 


54-127 


101 -204 


10.2 


16.0 


5.9 


3 


1.9-5.8 


70-120 


96-186 


8.8 


16.7 


14.7 


4 


1.3-4.3 


27-106 


103-198 




12.3 




5 


2.4-5.0 


65-118 


91 - 198 


10.2 


8.0 


9.7 


6 


2.3-5.6 


52-108 


109-177 


17.1 


1.4.7 


12.4 


TOTAL 


1.3-5.8 


27-127 


91-204 


14.0 


15.0 


10.3 



[0065] 

Table 2. Correlation coefficient (p) matrix. ABP is arterial blood pressure. See Table 1 
10 caption. 





CO 


MAP 


HR 


"RADIAL" ABP CO ERRORS 


0.11 


0.47 


-0.08 


FEMORAL ABP CO ERRORS 


0.09 


0.32 


0.14 



1 5 [0066] Since we were not able to obtain gold standard t d trends from the central ABP 
signals, we compared the two x D trends estimated from the "radial" and femoral ABP 
signals, which should, in theory, be equivalent. Table 1 summarizes the comparison results 
in terms of t d RMSNE in which the gold standard trends were established as the average of 
the two x D trend estimates. That is, the x D RMSNE here equals the root-mean-square of the 

20 difference between the two t d trend estimates divided by their sum and is the same for either 
signal. The low overall t D RMSNE in the table is buttressed by a high overall p value 
(mean±SD) between the two x D trend estimates (0.85±0.08). Thus, there is solid agreement 
between the two t D trend estimates despite substantial differences in short-time scale 
morphology between the "radial" and femoral ABP signals (Figure 6), and the overall CO 

25 measurement error is only 14.6%. Preliminary studies show that the six-minute intervals 
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used here may be reduced to, for example, three-minute intervals without materially 
affecting the results. 

[0067] In certain other embodiments of the invention it may be desirable to adjust the 
proportionality constant AC to vary depending on the arterial blood pressure or heart rate, if 
5 it is found that there is a systematic difference between dependence of the estimated cardiac 
output and the absolute cardiac output on heart rate or arterial blood pressure. 
[0068] The foregoing description is to be understood as being representative only and is 
not intended to be limiting. Alternative systems and techniques for implementing the 
methods of the invention will be apparent to one of skill in the art and are intended to be 
10 included within the accompanying claims. 
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